Directed excitation transfer in vibrating chains by external fields 
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We study the coherent dynamics of excitations on vibrating chains. By applying an external field and matching 
the field strength with the oscillation frequency of the chain it is possible to obtain an (average) transport of an 
initial Gaussian wave packet. We distinguish between a uniform oscillation of all nodes of the chain and the 
chain being in its lowest eigenmode. Both cases can lead to directed transport. 
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I. INTRODUCTION 

The transport of energy or charge is fundamental for a 
large variety of physical, chemical, and biological processes. 
One of the most prominent examples is the energy transfer in 
the light-harvesting complexes in photosynthesis Qj. There, 
the energy of the captured solar photons is transported via a 
molecular backbone to the reaction center where the energy is 
transformed into chemical energy. Recent experiments have 
shown that coherent features of the transport process might be 
crucial for a high efficiency Gill). Usually, the system and 
the dynamics of the excitation (exciton) is modeled by open 
quantum systems where the system of interest, e.g., the light- 
harvesting complex, is coupled to an external environment. It 
has been shown that the environment can also support the co- 
herent dynamics 0]-[9]. 

Most of the models assume a time-independent Hamilto- 
nian motivated by the fact that, indeed, the network of chro- 
mophores underlying the energy transfer is rather static, even 
at higher (room) temperatures. However, this need not be the 
case. One can easily imagine the situation where the underly- 
ing molecule is not static but performs some kind of mechani- 
cal oscillation. Asadian et ai. have shown that certain types of 
motions can enhance the transfer efficiency when compared 
to the static situation JTUJ. In a related model, Semiao et al. 
studied the modulation of the excitation energies of coupled 
quantum dots driven by a nanomechanical resonator mode, 
also enhancing the transport efficiency ifTTl . Vaziri and Plenio 
showed that the periodic modulation of ion channels leads to 
the emergence of resonances in their transport efficiency lfl2l . 

Another influence on the dynamics can be external fields. 
Hartmann et al. have shown for the coherent transport of an 
initial Gaussian wave packet on a discrete (static) chain of 
nodes that by suitably switching the direction of a constant 
external field, one can achieve directed transport 1 13 1. There, 
the switching frequency has been matched with the Bloch os- 
cillation frequency. The effect of Bloch oscillations on the 
trapping of excitations has been studied by Vlaming et a/., 
finding that the trapping efficiency crucially depends on the 
strength of the external field (the bias) lfl4l[T5l . 

Clearly, mechanical motions and external fields are not re- 
stricted to energy transfer in molecular aggregates. Other ex- 
amples include cold atoms in optical lattices whose spacings 
can be periodically modulated [ 16 ] or waveguide arrays where 
the "external field" is achieved by a linear variation of the ef- 



fective refractive index across the array, see, e.g., IfTTl . 

A question we address in this paper is whether it is possi- 
ble to engineer the excitation transport in systems performing 
mechanical oscillations with a constant external field such that 
also here one obtains directed transport. 



II. MODEL 

We consider the excitation dynamics on a finite chain of 
N nodes with time-dependent couplings J n (t) between two 
adjacent nodes of the chain. The Hamiltonian in the node- 
basis {|n), n = 1..N} reads 

N N-l 

H^ = Y J E n \n){n\+ £ J n (t)(\n){n + 1| + |n + l)(n|), 

n— 1 n— 1 

(i) 

where the E n are the site energies. Now, in addtion we apply 
an external field with strength /, such that the total Hamilto- 
nian for an excitation on a vibrating chain reads 

N 

H s =H + fY / n\n)(n\. (2) 

71=1 

For chains whose nodes (molecules or atoms) interact via 
dipole-dipole forces, the couplings decay with the third power 
of the distance between the nodes. For fairly large distances 
between adjacent nodes the coupling to next-nearest neigh- 
bors can be neglected such that the assumption of nearest 
neigbor couplings in Hq can be justified. 

Now, there are two competing effects: On the one hand 
excitations in a static chain with an external field perform 
Bloch oscillations flT3] [TSl |T9l . On the other hand the time- 
dependent couplings can cause an enhanced transport effi- 
ciency ifToUL?! . If the oscillations are periodic, the distances 
between two adjacent nodes vary in a given interval. Short 
distances means stronger couplings and thus faster transport 
from node to node. Longer distances lead to weaker cou- 
plings and slower transport. Therefore, matching the Bloch 
frequency with the frequency of the chain oscillation should 
lead to an effective transport in one direction along the chain. 
The reason is that in the first half of the Bloch period Tg the 
distances between the nodes are smaller while in the second 
half of Tb the distances are larger. This leads to different 
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displacements in the two half periods and consequently to an 

overall displacement of the inital excitation in one direction. 
We choose two scenarios for the couplings J n (t): 
(i) Each node of the chain oscillates uniformly with the 

same frequency u and with the the same amplitude a. The 

couplings follow now as 



Jn{t) = J{t) = -V/ 1 - 2asm(ojt 



(3) 



The same setting has been used by Asadian et aJ [ 10 1. In the 
following we also assume all site-energies to be the same, i.e., 
we set E n = E = 0. 

(ii) The chain is in its lowest eigenmode, such that for the 
gth eigenmode the couplings J nq (t) between nth and (n + 
l)st node are 



■/„,,(*) = -V/ 



1 — 2a n „ sin(ui a t + . 



(4) 



see Sec. lIIICl for details. 

Clearly, for time-constant J n — J we recover the known 
Bloch oscillations with frequency ujb = f /K (we set H = 1 
in the following). Thus, the period of the oscillation Tg = 
2tt/u b = 27!-//. 

The dynamics of an initial excitation is governed by the 
Liouville-von Neumann equation for the density operator 
p(t). Without any external environment leading to decoher- 
ence, the dynamics is fully coherent following 



A. Static chain 

We start by considering the static chain, i.e., no oscillations 
(a = 0). Without any external field and no external environ- 
ment, the dynamics of wave packets on the static chain is very 
similar to the motion of a quantum particle in a box 12T1I221 . 
One can also observe (partial) revivals of initially localized 
wave packets caused by reflections at the end of the chain, 
thus obtaining the discrete analog of so-called quantum car- 
pets J23] [24). 

When applying an external field, the situation changes. Fig- 
ure [T] shows for A = 78 the well-known Bloch oscillations 
in the occupation probabilities pkk{t) with Bloch frequency 
lub = f for / = 0.2 with no external coupling, A = 0, 
(left panel) and with small external coupling, A = 0.05, 
(right panel). One clearly recognizes the oscillation period 
of Tb = 2tv/ f = lOir. The coupling to the environment leads 
to a spreading of the wave packet over more and more nodes 
as time progresses. Eventually, this will lead to the equilibium 
distribution. 
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p(t) = -i[H s> p(t)}. 



(5) 



Now, if the system is coupled to an environment such that the 
total Hamiltonian can be split into three parts, H to t — H$ + 
Hr+ Hrs, where Hr is the Hamiltonian of the environment 
and Hps is the Hamiltonian of the system-environmental cou- 
pling. For small couplings to the environment we will study 
the dynamics by the Lindblad quantum master equation EU1 

N 

p(t) = -i[H s ,p(t)} (P(*) - <J'|p(*)|j»li>01. (6) 



where we assumed Lindblad operators of the form VX\j){j\. 
The term proportional to A mimicks the influence of the en- 
vironment leading to decoherence. In the following we will 
consider the occupation probabilities pkk(t) = (k\p(t)\k) for 
a given initial condition p(0). 



FIG. 1: (Color online) Static chain: Contour plot of the occupation 
probabilities for a = and / = 0.2 with A = (left panel) and 
A = 0.05 (right panel). Dark (black) regions correspond to large 
probabilities, while bright (yellow/white) regions correspond to low 
probabilities. 

In a continuous approximation for an infinite line, the po- 
sition of the center of the wave packet follows for vanishing 
initial momentum as lfl3l [T8l [T9l 



AN(t) = N c (t) — iVrj — — 



W_ 

T' 



! (/V2). 



(7) 



Obviously, there is no transport after integer values of Tb, 
only after Tb/2 = ixjf has the wave packet travelled by 
\AN(T B /2)\ = W/f = 20 nodes in the direction of the 
field. We note that by instantly reversing the field after Tb /2 
the wave packet will continue to move to the left side, such 
that it is possible to obtain directed transport by switching the 
field every half -period, see [13] for details. 



III. RESULTS 

In all calculations shown below we used N — 103 and an 
initial Gaussian wave packet centered at N C (Q) = N with 
a standard deviation of a = 6. We adjust Ao such that in 
the first two periods of the Bloch oscillations the wave packet 
does not encounter the edges of the chain, such that we can 
exclude interference effects caused by reflection. We further 
take V = 1. 



B. Uniformly oscillating chain 

If the chain is not static (a ^ 0) but oscillates such that the 
couplings are given by Eq. ([3]), it is possible to obtain - on 
average - a net transport of the wave packet in one direction. 
However, this will depend on the choice of the field strength 
/, i.e., on the frequency of the Bloch oscillation, on the phase 
shift 4>, and on the amplitude a. 



3 



1. Analytical approximation 

Before turning to the numerical results, we give an analyt- 
ical estimate of the displacements A A) = [N c (ITb) — No] 
{I E N) of the center of the wave packet after integer values 
of the Bloch period Tb- For the static (infinite) chain, starting 
from Eq. (|7]l and differentiating with respect to time, one has 

N(t) = -4Vsin(ft/2)cos(ft/2) = -2Vsin(/t), (8) 

which gives the temporal change of the displacement. Thus, 
the rate of transport from node to node is V. We extend this 
idea to the oscillating chain and replace the coupling V with 
the time-dependent coupling Jit). Then, we define the ap- 
proximate displacements by integrating N(t) over integer val- 
ues of Tb'. 



AN, 



(,approx 



-IV J dt 
o 



S m(ft) 



[l- 2a sin(ut + </>)] ' 



(9) 



For ui = f this leads to 



AA^/ a p plox 



YllixaV cos <fi 
f(l - 4a 2 ) 5 / 2 



(10) 



Clearly, the displacement is maximal for cf> = and minimal 
(zero) for <f> = tt. Note that AA) iapprox is only valid for the 
infinite chain. In the following we will compare AA^ approx to 
numerical results obtained from Eq. (|6j. As we will show, for 
the uniformly oscillating chain, AA). approx agrees very well 
with the numerical results. Also for the chain in its lowest 
eigenmode we will use AA^ approx as a starting point to define 
an adhoc fitting function AA^gt which also turns out to be in 
very good agreement with the numerical results. 



2. Numerical results 

Figure [2] shows the occupation probabilities p k t (t) for the 
case uib = f = oj = 0.2 with a = 0.1 and for different phase 
shifts qt>. Again, the left panels show the results for isolated 
chains (A = 0) and the right panels for small couplings to 
an external environment (A = 0.05). Plots in different rows 
correspond to different 0. Matching / with ui and having no 
phase shift results - on average - in a directed transport of the 
initial wave packet in the direction of the field. In the second 
half of each Bloch period Tb the wave packet moves in the 
opposite direction. However, this is overcompensated by the 
motion in the direction of the field in the first half of each 
period. 

The dependence on the phase shift can be expressed by only 
considering the average displacement AAV Figure [3] shows 
the dependence of AA^i and AN2/2 on <j> for the same pa- 
rameters as in Fig. [2] but with Nq = 52. Changing the initial 
condition to the center of the chain allows to vary in between 
and 2ir and thus avoiding interference effects due to reflec- 
tions at the ends of the chain. Note that this has no influence 
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FIG. 2: (Color online) Oscillating chain: Contour plot of the occupa- 
tion probabilities pkk{t) for a — 0.1 and w = / = 0.2 with A = 
(left panels) and A = 0.05 (right panels). The three rows correspond 
to different values of <f> = 0, n/4 and n/2, respectively. 



on the dynamics because the couplings in the chain are trans- 
lational invariant. We distinguish between AA^i after one and 
AA^ after two periods because, in general, one cannot expect 
a linear behavior of ANi in I. However, as it turns out ANi is 
approximately linear in I for the uniformly oscillating chain. 

Changing the phase shift allows to control the transport: No 
phase shift (0 = 0) results in values of AA^ w —21 after one 
period. A phase shift of <j> — n/2 results in a behavior similar 
to the Bloch oscillations in the static chain, i.e., no transport, 
see also Fig. [T] Increasing qt> further leads to a reversed mo- 
tion, i.e., the wave packet moves "uphill" against the direction 
of the field. For tj> = tt the maximal displacement after one 
period of AA^i sa 21 is obtained. For the uniformly oscillat- 
ing chain, the values for AN 2 /2 coincide with the ones for 
AA^i leading to the linear behavior A A) = lANi. In ad- 
dition, Fig. [3] shows the analytical estimate of Eq. (|9]l which 




FIG. 3: (Color online) Oscillating chain: Displacements ANi/l with 
/ = 1,2, extracted from the occupation probabilities for TVo = 52, 
as a function of (f> for a = 0.1, ui = f = 0.2 and A = 0. The dashed 
lines show AAT; ia p prox given in Eq. d9). 
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FIG. 4: (Color online) Oscillating chain: Displacements ANi/l with 
Z = 1,2 as a function of a for iV = 52, w = / = 0.2, = 0,7r/ 2 ,7r 
and A = 0. The dashed lines show AiV;, a p prc>x given in Eq. 



agrees with the numerical results. 

The magnitude of the displacements AJVj also depends on 
a. Figure p] shows AiVj/Z as a function of a for Nq = 52 and 
4> = 0, 7T/2, and 7T. While for <j) — n/2 there is no displace- 
ment after integer values of Tg, the displacements for = 
and for = 7r grow with increasing a. Again, the dashed lines 
show the approximation AiVi a p prox which nicely agrees with 
the numerical results. 

The effect of having directed transport depends on having 
the field strength in resonance with the chain oscillation fre- 
quency. In order to see how crucial the exact matching of / 
and uj is, we study slightly detuned frequencies u>b, i.e., a mis- 
match between uj and /. Figure [5] shows for <fi = (leading 
to maximal AiV; for uj = /) and for uj — 0.2 the occupa- 
tion probabilities pkk(t) for different values of /. Note that 
changing / also changes the Bloch period Tb = 2?r//, thus 
the time axes are different for different /. A field strength of 
/ = 0.18 or / = 0.22 reflects a detuning by ±10% of u. This 
still results in an average directed transport after two periods 
of AiV 2 w 30 for / = 0.18 and of A7V 2 « 36 for / = 0.22. 
Increasing the detuning further diminishes the transport. 

Figure [6] shows the displacements ANi/l for = and 
different values of /. The maximal displacement is obtained 
for uj sa /, as expected. Decreasing or increasing / results in 
smaller displacements: For f > uj the decrease in displace- 
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FIG. 6: (Color online) Oscillating chain: Displacements ANi/l with 
I = 1,2 as a function of / for a = 0.1, uj — 0.2, = 0, and 
A = (note the semi-logarithmic scale). The dashed lines show the 
approximation obtained by numerical integration, see text for details. 



ment is slower than for / < uj. One also observes that the 
displacements change directions. For / < uj, AA^/2 changes 
direction at about f /uj = 0.8 and ATVi at about f /uj m 0.67. 
For / > uj, the direction change happens at larger deviations 
from the resonace condition. Additionally, there are maximal 
displacements in the opposite direction. 

As before, we can obtain an approximation to the numerical 
results: Considering now / ^ uj in Eq. ^ and numerically 
integrating over integer values of the Bloch oscillation yields 
the dashed curves shown in Fig. [6] Again, the approximation 
is in very good agreement with the numerical data. 

Having now explored a large region of the parameters / /uj, 
a, and <\>, we see that the dynamics of an initial Gaussian wave 
packet can be manipulated by a suitable choice of these pa- 
rameters: We can make the wave packet move - on average - 
in one preferred direction by choosing the phase shift <fi. The 
magnitude of the displacements in either direction is given 
by a. Moreover, we do not have to exactly match the Bloch 
frequency ujb = f with the oscillation frequency uj in order 
to obtain directed transport, there is a fairly large range of 
roughly ±10% around f/oj = 1 in which large displacements 
can be obtained. 




FIG. 5: (Color online) Oscillating chain: Contour plot of the occupa- 
tion probabilities pkk(t) with uj — 0.2, a = 0.1, = 0, and A = 
for different values of /. 



C. Chain in lowest eigenmode 

In contrast to the previous section, we now consider the 
dynamics on a finite chain in its lowest eigenmode. Although 
this mode is similar to the uniform oscillation, the finite size of 
the chain becomes crucial leading to a non-uniform oscillation 
of the nodes. 

The couplings J n ,q{t) in Eq. Q between the nodes are ob- 
tained from a normal mode analysis of a free chain of nodes 
connected by springs, see [ 10] for details. Although the mo- 
tion of the nodes is not uniform E5ll . there are close similari- 
ties to the results presented in the previous section. 

In order to obtain comparable results we have to adjust the 
amplitudes and frequencies according to the couplings J n , q (t) 
between nodes n and n + 1 for the qth eigenmode. The cou- 
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plings in Eq. Q can be written as 

2asin[2u;£sin(g7r/2jV) + 0] 



-V 



cos(qn/2N) 



x sin(ng7r/iV) sin(<77r/2Af) 



such that one has 



and 



a sm(nqn/N) sm(qir/2N) 
cos(qw/2N) 



uj q = 2lj sm(qn/2N). 



(11) 



(12) 



(13) 



Thus, in the following we will use u) q = f as the resonance 
condition for the frequency and the field. For the amplitude 
a n>q to be comparable to the amplitudes in the previous sec- 
tion, we consider the average absolute value of the amplitudes, 
i.e., 



N 



N 



(14) 



TV 

n=l 

^ tan( g7 r/2iV) cot( g7 r/27V) = — 



aq 



Thus, we consider amplitudes a q which - on average - are of 
the same order as the ones in the previous section. This means 
that we choose the parameter a in Eq. ( 12 1 to be a = Na q /q. 

Similarly to Fig. [2} Fig. [7] shows the occupation probabil- 
ities pkk{t) for the case ub = f = wi. All plots in Fig. [7] 
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FIG. 7: (Color online) Lowest eigenmode: Contour plot of the occu- 
pation probabilities Pkk(t) f° r Si — 0.04 and Wi = / = 0.2 with 
A = (left panels) and A = 0.05 (right panels). The three rows 
correspond to different values of <j> — 0, 7r/4 and n/2, respectively. 
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FIG. 8: (Color online) Lowest eigenmode: Displacements ANi/l 
with I = 1, 2 for A r o = 52 (upper panel) and for No = 42 and 
No = 78 (lower panel) as a function of for ai = 0.04, Wi = / = 
0.2, and A = 0. The dashed lines show the fits for AAT;,Ht given by 
Eq. fB). 



show results for Si = 0.04. We use a% — 0.04 because 
this clearly avoids interference effects due to reflections at the 
ends of the chain. Coupling this system to an external envi- 
ronment leads, again, to decoherence and a spreading of the 
initial wave packet. 

Figure [8] shows a comparison of the displacements AJVj as 
a function of the phase shift for different Nq. Already for 
the central initial node, Nq = 52 (upper panel), one notices 
the asymmetry between the behavior of A N\ and AN2 /2 for 
values of G [0, 7r/2] and values of e [vr/2, tt]. For > it 
the difference between AJVi and AA^/2 is smaller than for 
< 7r/2, see in particular the points for = and <p = tt. 
One also notices that = 7r/2 yields ANi ^ 0, in contrast to 
the uniformly oscillating chain. However, the overall behav- 
iors for the two chains are very similar. Therefore, we fit our 
numerical result for A Ni by a cosine, as suggested by Eq. (|9j, 
namely, we use 



a at _ a aicos(0 + a;) 
AAT ijfit = ft_ 4 _ ?)5/2 , 



(15) 



where a; and /3; are (/-dependent) fit parameters. This already 
yields a very good agreement with the numerical results, see 
the dashed lines in Fig. [8] 

Changing the initial node Nq influences the behavior of the 
wave packet. Figure [8] also shows the behavior of AiVi and 
A A^/ 2 for N = 42 (lower panel, right half) and N = 78 
(lower panel, left half). While for N = 52 one has | AA^ | > 
|AAT 2 /2|, one observes for N = 42 and for iV = 78 that 
I A Ai| < \AN2/2\. However, for all initial nodes shown in 
Fig. [8] the maximal displacements (for = and = tt) are 
in the same region about \ANi/l\ m 12. 
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The slight asymmetry can be attributed to the non-uniform, 
i.e., non-translational invariant, motion of the nodes of the 
chain and the additional influence of the external field, which 
breaks the point symmetry with respect to the center. 
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FIG. 9: (Color online) Lowest eigenmode: Displacements ANi/l 
with I = 1, 2 as a function of Si for Wi = / = 0.2, and A = 0: 
upper panel for iVo = 52 with </> = 0, 7r/2, 7r and lower panel for 
No = 42 with 4> = tt/2, tt and iVo = 78 with = 0, vr/2. The 
dashed lines show the fits for AA^^t given by Eq. \15) . 

The Si-dependence of the displacements is shown in Fig. [9] 
Although the absolute values of ANi jl are different for dif- 
ferent No, there is a similar behavior for different values of 
(j>. Moreover, the behavior is similar to the one for the uni- 
formly oscillating chain, see Fig. |4] Therefore, we fit the in- 
dependence of ANi by AiVi fit given in Eq. ( 15 I. Also here 



are the fits in very good agreement with the numerical results. 

Figure 10 shows the displacements ANi and AN 2 /2 as a 
function of / for cf) = 0. Similar to the oscillating chain, the 
displacements are maximal for / k wj. The dashed lines 
show the approximations obtained for the oscillating chain 
(see Fig. [6|l but rescaled by a factor 1/2. Already this rough 
approximation yields good agreement to the numerical results. 
However, the points for f /to = 0.6 have to be considered with 
care, because such a detuning leads to interference effects due 
to reflection at the end node of the chain after one half period. 
This interference obviously can influence the dynamics of the 
wave packet. 

Now, also for the chain in its lowest eigenmode we obtain 
similar results to the ones for the oscillating chain. However, 
the absolute values of the parameters are different. Never- 



theless, the approximations given by Eq. (|9]l turn out to give 
qualitatively the correct behavior. Therefore, the same con- 
clusions as for the oscillating chain apply here. 
5r 




FIG. 10: (Color online) Lowest eigenmode: Displacements ANi/l 
with I = 1, 2 as a function of / for Si = 0.04, ui = 0.2, = 0, 
and A = (note the semi-logarithmic scale). The dashed lines show 
the approximations of Fig. [6] scaled by a factor of 1/2, see text for 
details. 



IV. CONCLUSIONS 

We have studied the coherent transport of excitations on a 
finite chain with time-dependent couplings between adjacent 
nodes of the chain and in the presence of an external field. The 
field leads to Bloch oscillations while regular time-dependent 
couplings can lead to an increased transport efficiency of exci- 
tations along the chain. We showed for uniformly oscillating 
chains and for a chain in its lowest eigenmode that match- 
ing the Bloch oscillation frequency with the frequency of the 
chain leads to an (average) directed displacement of an initial 
Gaussian wave packet. Applying a phase difference allows to 
manipulate the direction of the transport, while changing the 
amplitude of the regular oscillation allows to manipulate the 
strength of the displacements. We corroborate our findings by 
an analytic (continuous) approximation for the average dis- 
placement of an initial Gaussian wave packet in an infinite 
chain after integer values of the Bloch period. For the uni- 
formly oscillating chain, this ansatz yields a functional form 
for the displacements, which agrees very well with the nu- 
merical data. Using the same functional form also allows to 
define a fitting function for the chain in its lowest eigenmode, 
also leading to very good agreement with the numerical re- 
sults. In both cases, interference effects due to reflections at 
the ends of the chains have been neglected. 
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